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1. Introduction 

From the early development of modified (higher order) gravity models pQ, their 
dynamical equivalence to different representations has been an issue of research and 
intense debate. What the term "dynamical equivalence" means is that a particular 
gravitational action can be re-expressed as a new one, with a new set of field variables, 
and that there is an invertible mapping that relates the two sets of field variables, as 
well as the two actions (or Lagrangians) respectively. The variational principle of the 
new action will in principle require different boundary conditions and possibly different 
Gibbons-Hawking terms as well. Probably the most well known example of such an 
equivalence is that of the f(R) action to Brans-Dicke and Einstein-Hilbert one, the 
first through the introduction of an auxiliary scalar field [2] [I], an d the second through 
a conformal transformation of the metric [21 SI E] , [8] |12j. 

The motivation of re-expressing a gravitational action by introducing a new set of 
field variables might be related mainly to two things. The first is mathematical simplicity 
and convenience, if the new set of variables is to make the calculations one wants to 
perform simpler. Furthermore, if one is able to move from a variational principle that 
will lead to fourth order equations of motion to another that will lead to second order 
ones, that could be a benefit, since second order equations are in principle easier to 
handle as well as to interpret physically. The second possible motivation is related 
to physics. One gets more intuition and understanding of a gravitational theory, by 
studying its equivalence to other ones (like the equivalence between f(R) and Einstein 
gravity) . 

From a physical point of view care must be taken in the interpretation of physical 
quantities in the two different representations. The question that often arises is which 
of the two representations is the physical one. For example, in the case of the conformal 
equivalence between f(R) and Einstein gravity, the inclusion of matter in the action 
can raise the question of along which of the two metrics (original and conformally 
transformed one) do matter particles actually fall. For some interesting discussions on 
the subject one can refer to Refs [T3] -[2"3]. 

Given the equivalence between the bulk parts of two actions, this does not mean 
that the equivalence holds for the surface parts as well. More precisely, given the 
GH term of an action in one representation, then the GH term calculated using the 
equivalence with the other reresentation is possible to be problematic as far as the 
particular variational principle is concerned. As we will discuss later on, this is the 
case for f(R) and R + f{G) gravity, when the equivalence between the original and the 
Jordan frame action is considered. As we will see, the latter equivalence breaks on the 
boundary surface. 



In this paper, we will focus on f(R) [251 ESI EH EE] and R + f(G) [23 E8 ED] 
models, where G is the Gauss-Bonnet (GB) term defined as [31] 



"/" being an in principle non linear function of its arguments. We will study the 



G = i? 2 — AR a ^ R a p + R a p^ysRa^^ i 
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dynamical equivalence of above theories to other representations in vacuum, using as 
our main tool the rather general approach of Legendre transformation, which for the 
case of f(R) coincides with the standard procedure of introducing an auxiliary scalar 
field followed by a conformal transformation, something which is not true for R + f(G) 
theories, i.e. the latter cannot be conformally transformed to a minimally coupled, 
scalar-tensor frame. However, for the R + f(G) theory, we will show how using a 
Legendre transformation we can re-express it as a second order theory, with a new extra 
rank two tensor field. We will work at the level of the action and we will include in 
our analysis the relevant Gibbons-Hawking terms [32] . which are important for the 
consistency of the initial value formulation of the theory. Furthermore, we will calculate 
them explicitly wherever necessary. 

The structure of the paper is as follows: In Section |2] we briefly describe the 
approach of Legendre transforming a higher order gravitational action. In view of the 
latter approach, in Section [3] we show the equivalence of the full (including the relevant 
GH term) f(R) action to the Einstein-Hilbert one. We also discuss the equivalence 
between the relevant GH terms. Then, in Section [5] we consider the R + f(G) action 
and after calculating the GH term in the Jordan (scalar-tensor) frame, we explicitly 
study the effect of conformal transformation on the full Jordan frame action. Finally, 
in Section [6] we re-express the original R + f{G) action with one scalar and two extra 
tensor fields present (apart from the metric), one of them though being independent, 
and then discuss the classical dynamics of the system. Useful formulas and calculations 
can be found in the Appendix. We will work in natural units, c = G = 1. 

2. Establishing dynamical equivalence 

In this section we will briefly describe the idea and motivation behind Legendre 
transforming a gravitational action as a tool of moving to a new, dynamically equivalent 
variational principle of second-order, first applied in Refs [33| 131] . 

The higher order nature of non-linear gravitational Lagrangians like the f(R) or 
R + f(G) ones, comes from the fact that they are non-linear with respect to the second 
derivatives of the metric, since "/" is an in principle non-linear function with respect 
to its arguments, and R a pys — R a p-y5 {d 2 , (Vg) 2 , V 2 g). 

However, we can try to re-express higher order gravitational Lagrangians linearly 
with respect to V 2 g, by making them linear with respect to the curvature tensors |jj 
through the introduction of the appropriate "velocities" and "momenta" , in a similar 
fashion to the ordinary Hamiltonian formalism. The new variational principle will then 
lead to second order equations of motion for the new set of field variables. 

The appropriate identification for the generalised "position" and "moment" as well 
as the Legendre transformed Lagrangian will read schematically as 

q+*g*f), qt^ {R,R^,R a ^ s ), (2) 

| This is because curvature tensors are linear with respect to V 2 <?. 
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L(-g) 



1/2 



QiPi - Qi(q,p)Pi - (~g) 1/2 L(q,p) = (jiPi - H(q,p), 



(3) 



assuming the invertibility condition holds, d 2 L / \dq\dqj) 7^ 0. Quantities entering L 
will be in principle tensor objects. Different gravity actions give us different options in 
defining generalised "velocities" and "momenta". This will be made clear in Sections 
H] and El where we will apply the above formalism for the case of f(R) and R + f(G) 
gravity respectively. 

3. Dynamical equivalence of f(R) gravity, part I 

It is well known in the literature that through the introduction of an auxiliary scalar, 
the f(R) action can be re-expressed as a non-minimally coupled scalar-tensor one (also 
called Jordan frame action), and it is the latter that is usually conformally transformed 
to the so-called Einstein frame action. In this section, we we will focus and attempt 
to clarify the role of the relevant GH terms in the two representations, original and 
Jordan frame one. Then, in the next subsection we will demonstrate the equivalence 
of the full action (bulk and surface part) using the general approach of the Legendre 
transformation. Starting from the bulk f(R) action on a manifold M, 



J M 

through the introduction of an auxiliary scalar ip, one can re-express it in a dynamically 
equivalent way as 



with $ = f'(ip), and V($) = - / (J' -1 (<&)). For the latter we require that 

f"(ip) 7^ 0, so we are able to solve for ip = / /_1 (<I>). Action (jSJ) is the so-called Jordan 
frame action. 

The transition to the Einstein frame action will be shown explicitly in Section H] by 
means of a Legendre transformation. 

3.1. The f(R) Gibbons-Hawking term in the Jordan frame 

When considering gravitational actions on manifolds with boundary S, the variation 
gives boundary terms containing normal derivatives of the metric variation V« 5g a p- 
However, a well defined variational principle requires that only a particular set of 
dynamical coordinates (g a p and possibly its derivatives up to some order depending 
on the theory) is fixed on the boundary. In order to cancel the extra, unwanted surface 
terms, one needs to add a so-called Gibbons-Hawking (GH) term in the action [32]. The 
appropriate modification of the Einstein-Hilbert action turns out to be 




(4) 




(5) 




(6) 
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where h is the induced metric on the surface E, and K is the trace of its extrinsic 
curvature. Variation of (Q is then performed keeping only the metric g a @ fixed on the 
boundary. 

In the original f(R) action (J4]) there is no natural GH term which cancels the extra 
unwanted higher derivative boundary terms, except for the particular case of maximally 
symmetric spacetimes [35J. However, in the Jordan frame of f(R) one can find an 
appropriate Gibbons-Hawking term in full generality, and then under the assumption 
that the dynamical equivalence between different representations holds on the boundary 
surface, one is able to re-express it in the original f(R) representation. Let us examine 
this more carefully. 

First we want to find the GH term in the Jordan frame, and so we vary (jSJ) with 
respect to the bulk metric and after discarding the bulk contributions we get 

SS^ = - [ d n - l x V^hV p ($ g aP 6r p a p) = 2 I d n ~ l x V^h$5K. (7) 
Jt, Jt, 

Therefore, the GH term that should be added in the Jordan frame action ([9]) is 

S^ = -2 J d n ~ x x yf^h^K, (8) 

and variation should be performed with (Sg a p, 5<E>) vanishing on S. After using the 
correspondence between Jordan and original frame, $ <R- f'{R), we find the GH term 
in the original frame to be [§| 

S^ = -2 J d n ~ l x V Z hf'{R)K. (9) 

However, variation of this boundary term generates a new term 

<5S E = -2 J d n ~ x x V^hf"(R)K5R, (10) 

which vanishes only by requiring that 5R = on the boundary. In Ref. [36] it has 
been shown that the GH term ([9]) is necessary in order to derive the correct Wald 
entropy for f(R) gravity. However, keeping R fixed on the boundary surface can be 
problematic. Since R includes both the first and the second derivatives of the metric, 
keeping it fixed would in principle require that the second derivatives of the metric are 
held fixed too, which overconstrains the actual formulation. The only possibility that 
would prevent the latter from happening would be that the condition 5R = is satisfied 
through some special configuration of the field variations on the boundary, something 
that would restrict the generality of our variational principle. 

Therefore, we see that equivalence between the two representations breaks down 
at the boundary, when the consistency of the variational principle is considered. This 
failure indicates that the two theories cannot be truly equivalent. Furthermore, the two 
theories are inequivalent at the quantum level as well; considering the path integral 

§ We obtained the GH term in the original f(R) representation by substituting the equation of motion 
for <I>, V'($>) = R(g), into Eq. ((HJ. This makes clear that the equivalence is demonstrably valid only 
on-shell (i.e. at the level of the classical equations of motion). 
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defined in the Jordan frame, the integration over $, will generate extra terms in 
the effective action, making the latter inequivalent to the one defined in the original 
representation. 



4. Dynamical equivalence of f(R) gravity, part II 

We now want to exploit the dynamical equivalence between the full f(R) and Einstein- 
Hilbert action, using solely the Legendre transformation approach, presented in section 
El which is more general and for the f(R) case gives the same result with the conformal 
transformation. This was done in Ref. [33] for the bulk Lagrangian, discarding a total 
derivative. In the following we will show how to cure this by Legendre transforming the 
GH term, apart from the bulk part, getting this way the correct Einstein-Hilbert GH 
term as well, [jj 

Let us begin with the action 

jn—l 



S= / d n x^f{R) + / d n ~ 1 xV-hf'{R)K, (11) 
Jm Jt. 

including the GH term found in the previous section^. As we will see below, the inclusion 

of the latter is indeed a good choice. We will need separate variables for the boundary 

surface, and the Legendre transfomed action will be of the form 

jn— 1, 



S= d n x y/=jj (q B p B -H B )+ / d n ~ L x V-h - H s ) , (12) 

Jm Jt, 

with Hi = Hi(q,p). Let us first naively associate for the generalised bulk velocity, 

q -H- R a p. Then we get the bulk conjugate momentum as 

1 BT 

p B (q, q)^~g a ^^j^ = f\R)9 aP . (13) 

We see that the definition of the conjugate momentum defines a conformal relation 
between two different metrics. In fact, as we will see below, g a p is the metric in the 
Einstein frame. However, the correct association for q is not exactly R a p, but R, since 
relation ( IT3l) cannot be inverted for R a p. 
We proceed by identifying 

q B <+R, g E O K, (14) 
p B (q, q) $ = = f'(R), (15) 



-gdR(g) 

pe(?, q) $ = ~^Zt^ = f'(R). (16) 



-h dK{h) 

|| The transition to the Einstein frame by means of a conformal transformation, including the relevant 
GH terms, has been studied in Ref. [36 . 

As it was previously discussed, the GH term @ can be in general problematic, however it turns to 
be neccessary here in order to cancel extra surface terms after Legendre transforming the action. 
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Now, using the intuition gained from (1131) . and using (jl5j) - (jl6j) . we define the following 
relations for the bulk and surface metric respectively 

g"f> = $V=n} gO0 = $(^7/^. (17) 

The invertibility condition is not satisfied on the boundary £, since d 2 LY,/dK 2 = 0. 
However, the surface part of L can be still defined, with the only difference that the 
surface Hamiltonian will vanish identically, H% = 0. The bulk Hamiltonian is calculated 
after solving one of relations (US]) for R, H B (&) = Qf'^ 1 ^) - /(/ /_1 ($)). Using this 
together with (| 14ft - (jib]) , and substituting in ( lT2"j) we arrive at 



S 



(18) 



(Txy/=g[$R(g) - H B ($)] + / d^xV^h^K. 

*M JT, 

For the transition to the Einstein frame we will use the general equations (lA.lOj) 
and (lA.lip relating two Ricci (extrinsic curvature) tensors, evaluated for two different 
metrics g a p (h a/3 ) and g a p {h a p). Defining $ = exp[0] and uj{n) = (n — l)/(n — 2) we 
get 



5 



B 



d n x 



'M 

+ 2u{n) 



R{g) - u{n)d K ct)d K ct) - e(^)^(0) 



ft 



d 



n—l 



X 



-9 



(f-'xV -h(d K (j>)n h 



2K(h)-2u(n){d K 



(19) 



(20) 



and after summing up we arrive at 



S 



d n x 



M 



■ — ■ n i 

R(g) - u(n)d K <i)d K (t) - e^-^> H B ((f>) 



+ 2 d 



in— 1, 



X 



hK(h), (21) 



h K denoting the normal vector to E. We see that we arrive at the correct, full Einstein- 
Hilbert action, following a conceptually different and more fundamental procedure. 
The conformal relation between the two metrics was revealed naturally through the 
definitions of the conjugate momenta. 



5. Dynamical equivalence of f(G) gravity, part I 

In this section we will aim to express the Jordan (scalar-tensor) frame of the R + f(G) 
action as a minimally-coupled theory by means of a conformal transformation. Firstly 
we will derive the GH term in the Jordan frame, and then find the appropriate one in 
the original frame, as dictated by the equivalence between frames. Then, in subsection 
15.21 we will continue with conformally transforming the full, Jordan R + f(G) action. 
Our starting point is the R + f(G) action 

S= [ d n Xv ^[aR+f(G)], (22) 

with G defined in ([1]) and a a dimensionless constant. 
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Through the introduction of an auxiliary scalar field ip we get the Jordan frame 
action as 

S J = f d n x^[aR + $G - V{$)], (23) 
with $ = f(^), and V($) = - / assuming that f'(^) 7^ 0. 



5.1. T/ie /(G) Gibbons-Hawking term in the Jordan frame 

The motivation of this subsection is the same as in the f(R) explained in Section 

13.11 We will derive the appropriate GH term in the original action fl22|) as dictated by the 
equivalence with the Jordan frame, by first calculating the Jordan frame one, presenting 
the explicit results of the surface parts of the action variation. Some useful variation 



formulas and definitions used can be found in |Appendix A 
We start from the Jordan frame action 



and vary each of the Gauss-Bonnet 
terms separately with respect to g a p using relations flA.lj) - (I A . 6 1) . We focus on the f(G) 
term, since the GH term for R is given by ([H} for f'(R) = 1. We will again present only 
the boundary part of the variation, as well as work in Riemann and Gaussian normal 
coordinates [37]. With the aid of integration by parts, and using equation (1A.1I) . we get 



SS^ = -4 / d n ~ l x V-h $R5K, 



(24) 



ss 3 ^ 



I I d n - l x^h<$> 
- n a n p n K R KX V) 



2n p R aK V K - n x R a/3 V x - n x h a/3 R KX V K , 
8g a 8, (25) 



d 



n—l , 



X 



-h^[n x R a ^V K ] 5g aP . 



(26) 



The geometric relevance of the above terms becomes evident if we express them in terms 
of tensor objects defined on the boundary surface using the Gauss-Codacci equations 
Doing this, and adding up all three terms together, we arrive at 



5S- 



d 



n—l 



2KK^5K M + K 2 5K - K au K afi SK) 



(27) 



with G^ 1 the Einstein tensor defined on S. Since we require that 8g a p = on £ (or 
o~h a/ 3 = 0), it follows that 5G a p = and 5K a p = 5K alB = 5K a p on S as well. Using those 
facts, we can go backwards in f l27|) and check that it is the variation of the following 
quantity 



d 



n—l, 



2 ( 2G af} K aP + J 



(28) 



with J = ?^K P k K kX K\ p — KK k \K kX + \K^. The appropriate supplement for the initial 
scalar-tensor action is therefore equation ( 12"5|) with a minus sign instead. The GH term 
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for a simple Gauss-Bonnet action (L oc y/—gG) has been derived under more general 
assumptions in a braneworld context in Ref. [38], as well as in Ref. [39] using the 
calculus of differential forms. 

Now, as in the f(R) case, we can use the equivalence f'{G) $, to find the GH 
term in the original f(G) frame if the equivalence is to hold on the boundary, 



d^xV-hfXG) 



2 ( 2G a ^K a p + J 



(29) 



Now, variation of the action requires 8g a p = and SG = on E. The latter condition 
can yield a problematic variational principle for the same reasons discussed in Section 
13.11 Therefore, for R+ f{G) theories as well, true equivalence on the boundary surface 
is broken. 



5.2. Conformal transformation of the Jordan frame action 



We now want to study if the non-minimally coupled, full Jordan frame action ([29]) . can 
be decoupled from $ and written in a G+ scalar field form, similar to the f(R) case, 
using a conformal transformation of the metric. 

Before we start with the calculations, let us introduce a notation that will make our 
equations look shorter. So, only for the rest of this section, we shall define: Hi = (n — i) 
and r» = where n is the spacetime dimensionality, rij is not to be confused with 
the surface normal n a . 

We shall begin with the bulk term. Using the transformation formula ( 1A.8|) together 
with the conformal factor identification 

$V(»-4) = ^ (30) 



n 



and omitting the potential which transforms trivially, the action ( 1231) after the 
redefinition <3> = exp[0] becomes 



d n x 



M 



-g$G^ / d n x 

J M 

- 8r 4 n 3 [0 Ta/3 



g{G 



■] R a j3 - 2r 4 n 3 [3r 4 n 4 ^ 



+ 4r A n 3 n 2 



2n 5 n 3 + 



> a <f>>P 
nni 



+ 4r 4 n 5 n 3 n 2 



3ar^ 5 me (1 - r4n2) 

c K + 2ar A me {1 - r4n2 



t - 200 - ae (1 " r4n2) 



15 



+ Ar 4 n 3 n 2 00 - r 4 (2n 5 + tii 



□0 



- 4r 4 n 3 n 2 ( 



Ta/3 



(31) 



with □ = V K V K . 

Identification (130]) breaks down for n = 4, and in fact it is valid only for n > 5. This 
means that we are unable to decouple the scalar $ from the Gauss-Bonnet term unless 
n > 5. For n > 5, the GB term is minimally coupled to the scalar $ = exp[0], but 
there are new couplings between the derivatives of 0, the Ricci tensor and Ricci scalar. 
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In this case, action ( 1ST]) , plus the scalar potential term of ( I2"3"j) . describes a fourth order, 
non minimally coupled scalar-tensor theory. 

Let us now turn attention to the conformal transformation of the relevant GH term, 



calculated in Appendix B One can see that the variational principle requires that we 
impose apart from 5g = and 5<p = 0, the extra conditions V5g = and V<5</> = on E. 
The R + /(G) action cannot be expressed as a second order, minimally-coupled scalar 
tensor one, in contrast with f(R) gravity. 

6. Dynamical equivalence of f(G) gravity, part II 

The richer structure of the R + f(G) action gives us more options in identifying 
generalised velocities, compared to the f(R) one. In this section we want to take 
advantage of the latter fact, and re-express the original R + f(G) action as a new one 
with not only a new scalar, but with new tensor fields as well, by means of a Legendre 
transformation. The new variational principle will be of second order. 
Our starting point is the action 



S= f (Pxy/=g[aR+ f(G)]. 



(32) 



We proceed with defining our conjugate momenta as 

1 BT 

p x ^t> = = « + 2iV'(G), (33) 

\/—g oR 

1 BT 

P2 » ## = = -%f(G)K#> (34) 

y/—g OK a/3 

1 BT 

p 3 o vj>« = = 2f(G)R a W. (35) 

yj-g dR a j3 j5 

Defining $ = $(G) = f'(G), the inverse of the above relations read 

R = ^ (* - a) , (36) 
= (37) 

Rj hS (*,g,g 1 <r) = ^<r a /h ', (38) 

with g K K = g KX g K \ and f'(G) ^ 0. In fact, we will use g a p to raise and lower indices for 
the rest of the section. 

For the calculation of the Hamiltonian we will need to express the Gauss-Bonnet 
term in terms of the new fields (\1/, g a/3 , a a ^ s ). Using the inverse relations ( l36l) -( )38l) we 

get 

1 

with the function V defined as 

r(¥, g, a) = (* - a f - + tr/V^. (40) 



G = R 2 - AR^R^ + R^paR^ = 7^ r (^> ~9, <r), (39) 



The dynamical equivalence of modified gravity revisited 11 

Furthermore, we assume that we can invert relation (1391) El and express the Gauss- 
Bonnet term in terms of the function T as 

G = G'\T) = J(r), (41) 

so that 

f'(G) = f'(J(T))^F(T). (42) 
Using all the above, we can now calculate the Hamiltonian as 

ff(r(*, g, a)) = *J?(<I>, g, a) + g" f ', 9, ") + aJ^ s R"g lS (9, g, a) - (-g)- 1 ' 2 ^, g, a) 

= ^T)-/Wn). (43) 

Notice that the fields (\l/, g a/B , o~ a l3 ' yS ) enter implicitly in the Hamiltonian through the 
function T. 

The Legendre transformed action then reads 



S[V,g,g,o] 



d n x ^J^g 



M 



*R(g) + ~g aP R^{g) + o- a ^ 5 R a ^ s (g) - H(T(^,g,a)) 



(44) 



To get the equations of motion we vary the action S with respect to the four fields 

(y,g a e,a a W,g a p) to get 



jL = R(g)-2H t (V-a) = 0, 



(45) 



SS 1 
jz^p = Rap(g) + -H'g a p = 0, 



5S 



R a /3-ys(g) — 2H'a a ^s = 0, 



(46) 
(47) 



5S 



8g 



a;3 



l -V p V p ~g aP + \g aP V K V x r X ~ 2V K V x a K ^ x 



-2** 

- l -H> 
2 



~g KX Rnx{g) + o- K ^ v R\, v {g) - H{Y) 



^V>, XK - g Kia gV 



0. 



(4* 



with G al3 = R af} - \g ap R, H' = H'(T) = dH/dY and covariant derivatives V Q defined 
with respect to g a p. 

Variation with respect to g a p yields surface terms oc Vg a p. We want to keep only 
the fields fixed on S and not their derivatives, so we have to add in action (]44j) the 
following GH term 

§ 2 = 2^ d n ^x V^h (VK + ~g^K[,n P] + a a ^ 5 r« h n s] ) , (49) 
+ The necessary condition is that [(/' (G)) 2 G}' ^ 0, implying that /(G) ^ CiVG + G 2 . 
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with [A,B] = | (AB - BA). 

If we now contract equation (j4T|) with g Q 7 and add it to f|47|) we get the relation 

0" 7 a 7 /3 = = ~ T<7ay8- (50) 

The latter implies that CTq,^ can be expressed as some combination of g a p and plus 
some traceless part, while the trace of that expression should give (I50|) . To find the 
latter expression we can expand the Riemann tensor in terms of he Ricci tensor and 
scalar according to 

Ca/3-ys - a n (g a [sR~,]p + gpfrRsia) - b n Rg a [ 7 gs]i3, (51) 

with a n = , b n = ( w _ 1 - ) 2 ( n _ 2 ) an d C a ^s the Weyl tensor which is traceless in all its 
indices. After use of equation relation (I5T]) can be expressed as 

Ra{3~{6 = CafoS + (9a{S9i\P + 9/3[-Y9S\a) + ~^H' 'gg a [j95]l3, (52) 

and plugging the latter into equation f|47|) to substitute for the Riemann tensor, we get 
a relation between a^s, g a p and g a p 

o~ a p~i& = -^pC a i3js{g) + (g a [8g~{]f5 + g^gs^) + gga^gsip- (53) 

Combining equations ( )45|) and ( l46l) we can find a similar relation for \l/ 



= (54) 

with g = g a Pg a /3. One would like to be able to solve equation fl53|) for cr Q/ 3 7 5 = o- aj3l s(g, g). 
However, this is not in principle possible unless H' = constant (corresponding to the 
trivial case of f{G) = G) or Cap-^g) = 0. The latter case includes the case of the FRW 
spacetime or maximally symmetric spacetimes like the Minkowski one. In that case, all 
fields can be expressed in terms of g a p and g a p and we can get a solution for the latter 
ones by solving the appropriate system of second order differential equations, which we 
derive below. 

Now, we want to derive a system of evolution equations for the set of fields (g a /3, gap)- 
The first equation we will use results from equation ( 1461) after taking its trace once, 
together with some simple algebra. To get the second equation, we use relations (153"]) 
and (1541) together with the Ca^sig) — ansatz to express the last of the equations of 
motion, equation ( )48|) . in terms of g a p and g a p alone. This way we arrive at the new 
system of second order equations for the set of fields (g a p,g a p) 

G a fi = ~2 H ' (da? - ^g a ^g K Kj , (55) 



H' 



r n ~g K{a fK + s n gg<# + 2ag afi 



(56) 
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with H = H(T(g,g)) and the operator P = P{g) as well as the constants p n , q n , r n , s n 
defined in Appendix C To arrive at equation (1561) we have used the following relations 



1 

16 



4(n - 3) 



{n-2f 



n 



n — 3 



n 



(57) 



and 



a 



[n 



g afi ~g m ~g PP + (n 



4)g<> a f P 



n 



n 



~g~g aP 



(5* 



as well as o- a ^s(g,g) given by (153"]) with C a p 1 s{g) = 0. 

A look at the first equation of the new system, equation (|55p . shows that at the 
level of the equations of motion we can express the dynamics as GR, minimally-coupled 
to an effective energy-momentum tensor (the r.h.s of the equation) described by the 
spin two field g a p. 

There is one extra constraint the fields satisfy, that is the Bianchi identities. Since 
the l.h.s of equation ( )55|) is covariantly conserved, as dictated by the Bianchi identities, 
then the r.h.s should be as well, 



H ' ( g<*p - 2&tf0 



o. 



(59) 

have to satisfy, together 



The latter equation is a condition the set of fields (g a p,g a p) 
with the equations of motion. 

We will not seek solutions of the system described by equations fl55|) - fl56|) in this 
paper, leaving this for a possible future work. However, it is easy to see that Minkowski 
space is a solution for (g a p,g a p) = (f/a/3,0). 

We see that at the classical level of this representation there are two independent 
fields, (g a p,g a p), satisfying a system of second order equations together with a second 
order condition, the Bianchi identity. These equations should be classically equivalent 
to the original fourth order ones, as g a p is related to the second derivatives of g a p. It 
should be noted that there is no reason to expect complete equivalence of the quantum 
equations, as the measure of the path integral in the different representations can 
introduce new terms. 



7. Conclusions 



Legendre transformations are a fundamental tool to study the dynamical equivalence 
between different modified gravity actions, with the aim of understanding better the 
nature of the theories under study. When working the level of the action, a consistent 
analysis should take into account the appropriate Gibbons-Hawking (GH) terms (full 
action). Although in a general context there are no natural GH terms for both f(R) and 
R + f(G) actions, however one can define them considering the dynamical equivalence 
between two different representations of the particular action on the boundary surface, 
as it was done in Section 15. 1\ when we considered the equivalence between the original 
action and the Jordan frame one. However, the GH terms found through this procedure 
turn out to render the variational principle inconsistent. This is due to the fact that the 
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two representations are not equivalent at the quantum level, as pointed out in Section 

Due to the structural simplicity of the full f(R) action, the Legendre transformation 
yields in this case the same result as a conformal transformation of the original action. 
On the other hand, the R+ f(G) Jordan (non-minimally coupled) frame action cannot 
be re-expressed as a second order theory through a conformal transformation, despite 
the fact that the auxiliary scalar decouples from the Gauss-Bonnet term for dim > 5. 
The resulting theory is still of fourth-order, as was calculated explicitly for the full 
action in Section [5j 

However, the more complex structure of the R + f(G) action, allows one to re- 
express it, by means of a Legendre transformation, as a second order theory with extra 
tensor fields apart from scalars. In the new representation, it turns out that only two 
fields are the independent ones, the metric g a p and the rank two field g a p. At the level 
of the equations of motion, we are able to recover GR, sourced by an effective energy- 
momentum tensor, which is a function of g a p. Although the two representations are 
classically equivalent in vacuum without boundary, at the quantum level they differ, as 
integrating out the extra fields generates new terms in the effective action. 
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Appendix A. Definitions and formulas 

The n-dimensional bulk metric g a p induces an (n — l)-dimensional metric h a p on the 
boundary surface £ as 

h a /3 = 9a/3 ±n a np, (A.l) 

for a spacelike (+) and timelike (— ) surface £ respectively and n a is the normal to the 
surface £. h a p can be used as a projection operator from the tangent space to the bulk 
M to the tangent space to the boundary E at a point P . Particularly, for its action on 
n a it is h a ^np = 0. 

The extrinsic curvature K a p is an (n — l)-dimensional tensor that measures the 
"bending" of £ in the bulk spacetime M, and is defined as 

K aj3 = ^£ n h a p = V a £g = liW^/j = /i a 7 V 7 ra /3 , (A. 2) 

where "<£"' is the Lie derivative, a unit tangent to the geodesic congruences orthogonal 
to E, n 13 any other normal to E, and V a defined with respect to the bulk metric g a @. 
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The following variation formulas are used to calculate the GH terms presented in 
Sections [3j |5] and EJ They are evaluated using the special coordinate systems of Gauss 
and Riemann normal coordinates respectively [37]. We have 

1 



6T C 



;9 ap (Sg. 



v\p,n] 







SR a p lS = g aK (6g K [sw]p - o~gpi r> 6} K ) , 



8K = 5K C 



-n"h^V p 5 gcn , 



(A.3) 
(A.4) 
(A.5) 

(A.6) 



with [A, B] = | (AB — BA) . The variation of the Ricci tensor and scalar can be found 
beginning from (1A.4|) and calculating the variation of the appropriate contractions, for 
example, 5Rp$ = 5(g 1 a R a p^). 

If g a p >-> 9ap = Q 2 9ap we have [HI E] 



r = q 2 



R + 2{n- - n{n - l)Q- 2 g a ^Q ta Q v 



(A.7) 



G = fi 4 G - An&- 1 (tiRapti * - RU£l\ 



n- 1 



K a0 - tt 1 h a0 Q, K n K 



(A.8) 
(A.9) 



where in ( 1A.8I) we use the convention ni = (n — i), with n the spacetime dimension. 

For two different metrics g a p and g a p, defined on the same manifold M, and not 
necessarily conformally related, we have 



Ra/3-y (9) — Rafi-y ((j) ~ V^C al — V a C p^ + C K aJ C p K — C K p 7 C 
~ ~ 1 ~ 1 

K a p{h) — K a p(h) = -£ n h a p — -£ n h a p = —h^C^-ypn^, 



(A.10) 



(A.ll) 



with C a p 1 = \g aa (ypg~ ja + V^gp a — V a gp^). Particularly in Section HJ where g c 
&g a p = $), equations ( 1A.10I) and ( 1A.11I) are used in the form below 



R(g) = $ 



Rid) - j{n - 2)(n - 1)$~ 2 <9 K $<9 K $ + (n - ^V^ -1 ^) 



(A.12) 



K = h afi K a p = $ 1/2 
with n a = § x l 2 n a and K = h al3 V a hp 



K - -{n - l)^ -1 ^ 



(A.13) 
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Appendix B. Conformal transformation of the Gauss— Bonnet GH term 



Here we will present the conformal transformation of the Gauss-Bonnet GH term, (128 
For the two terms of f l28|) we get respectively 



d n ~ x x 
s 

+ n-£T x 
1 



h fi 4_n J 



K - K aB K a > 

3 /r> _cck\3 



{Q K n K ) - n 3 n 2 K Q~ 2 (Q K n K ) 



(B.l) 



d n ~ l x V^h <l>G al3 K a8 i-> / d^xV-h fi 4 ~ n $<! G K, 



+ n 3 tt~ 2 



2 



"2 



(B.2) 



Adding up terms (1B.1I) and (IB. 2 1) we get the GB GH term in the conformally 
transformed frame. However, boundary terms resulting by variation of action (l3Tj) with 
respect to g a/ 3 and $, will not be able to cancel with the GH subterms in ( IB. II) and 
(IB.2j) . as was the case in f(R). Consequently, we are left with terms proportional to 
first and second order derivatives of both the metric and scalar field on the boundary 
surface E, which should be held fixed in the initial value formulation, together with g a p 
and themselves, in order for the GH term to be zero in the total variation. 



Appendix C. Complementary definitions for Section EJ] 



Here we define the operator and the constants used in equation (1561) . The operator is 
defined as 



P. 



h'X 
(lis 



(a/3) 



'I - a % 



while the constants a n , b n , c n , p n , q n , r n , s n respectively 



n-2 



(n - l)(n - 2) 



_ 1 -2b n 
i c n — , 
1 - a n 



a 2 1 
Pn = a n - ^-(1 + b n ), q n = -[1 - a n (a n - 2b n )}, 

r n = 2[1 - ^{n - 4)], s n = ^[2a n (a n - 2b n ) - 1]. 



(CI) 

(C.2) 
(C.3) 
(C.4) 
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